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Abstract 

We derive a weak turbulence formalism for incompressible MHD. Three-wave in- 
teractions lead to a system of kinetic equations for the spectral densities of energy 
and helicity. The kinetic equations conserve energy in all wavevector planes normal to 
the applied magnetic field -Boe||. Numerically and analytically, we find energy spectra 
~ k"^ , such that n+ + n_ = —4, where are the spectra of the Elsasser vari- 
ables = V lb b in the two-dimensional case (/cy = 0). The constants of the spectra 
are computed exactly and found to depend on the amount of correlation between the 
velocity and the magnetic field. Comparison with several numerical simulations and 
models is also made. 
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1 Introduction and General Discussion 



Magnctohydrodynamic (MHD) turbulence plays an important role in many astropliysical 
situations (Parker 1994), ranging from the solar wind (Marsch and Tu 1994), to the Sun 
(Priest 1982), the interstellar medium (Heiles et al. 1993) and beyond (Zweibel and Heiles 
1997), as well as in laboratory devices such as tokamaks (see e.g. Wild et al. 1981; Taylor 
1986; Gekelman and Pfister 1988; Taylor 1993). A very instrumental step in recognizing some 
of the features that distinguished MHD turbulence from hydrodynamic turbulence was taken 
independently in the early sixties by Iroshnikov (1963) and Kraichnan (1965) (hereafter IK). 
They argued that the destruction of phase coherence by Alfven waves traveling in opposite 
directions along local large eddy magnetic fields introduces a new time scale and a slowing 
down of energy transfer to small scales. They pictured the scattering process as being 
principally due to three wave interactions. Assuming 3D isotropy, dimensional analysis then 
leads to the prediction of a k~^^^ Kolmogorov finite energy flux spectrum. 

However, it is clear, and it has been a concern to Kraichnan and others throughout 
the years, that the assumption of local three dimensional isotropy is troublesome. Indeed 
numerical simulations and experimental measurements both indicate that the presence of 
strong magnetic fields make MHD turbulence strongly anisotropic. Anisotropy is manifested 
in a two dimensionalization of the turbulence spectrum in a plane transverse to the locally 
dominant magnetic field and in inhibiting spectral energy transfer along the direction parallel 
to the field (Montgomery and Turner 1981; Montgomery and Matthaeus 1995; Matthaeus 
et al. 1996; Kinney and McWilliams 1998). Replacing the 3D isotropy assumption by a 2D 
one, and retaining the rest of the IK picture, leads to the dimensional analysis prediction 
of a kj^ spectrum (Bq = -Boe||, the applied magnetic field, A;|| = k • ey, kj_ = k — k\\ ey, 
k± = |k_L|) (Goldreich and Sridhar 1997; Ng and Bhattacharjee 1997). 

A major controversy in the debate of the universal features of MHD turbulence was in- 
troduced by Sridhar and Goldreich (1994) (hereafter SG). Following IK, they assumed that 
the small-scale MHD turbulence can be described as a large ensemble of weakly interacting 
Alfven waves within the framework of the weak turbulence theory. However, SG challenged 
that part of IK thinking which viewed Alfven wave scattering as a three wave interaction 
process, an assumption implicit in the IK derivation of the k~^^^ spectrum. SG argue that, 
in the inertial range where amplitudes are small, significant energy exchange between Alfven 
waves can only occur for resonant three wave interactions. Moreover, their argument con- 
tinues, because one of the fluctuations in such a resonant triad has zero Alfven frequency, 
the three wave coupling is empty. They conclude therefore that the long time dynamics of 
weak MHD fields are determined by four wave resonant interactions. 

This conclusion is false. In this paper, we will show that resonant three wave interactions 
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are non empty (see also Montgomery and Matthaeus 1995; Ng and Bhattacharjee 1996) 
and lead to a relaxation to universal behavior and significant spectral energy redistribution. 
Moreover, weak turbulence theory provides a set of closed kinetic equations for the long 
time evolution of the eight power spectra (to be defined below, in equations (21) and (22)), 
corresponding to total energy e*(k), poloidal energy $''(k), magnetic and pseudo magnetic 
helicities R^(k.), ^^0^) constructed from the Elsasser fields = v + sb, s — ±1, where v 
and b are the fluctuating velocity and Alfven velocity respectively. The latter is defined 
such that b = B/ ^/^oPo, where po is the uniform density and /^o the magnetic permeability. 
We will also show that a unique feature of Alfven wave weak turbulence is the existence of 
additional conservation laws. One of the most important is the conservation of energy on all 
wavevector planes perpendicular to the applied field Bq. There is no energy transfer between 
planes. This extra symmetry means that relaxation to universal behavior only takes place as 
function of k± so that, in the inertial range (or window of transparency), e*(k) — f{k\\)kj^ 
where f{k\\) is non universal. 

Because weak turbulence theory for Alfven waves is not straightforward and because of 
the controversy raised by SG, it is important to discuss carefully and understand clearly 
some of the key ideas before outlining the main results. We therefore begin by giving an 
overview of the theory for the statistical initial value problem for weakly nonlinear MHD 
fields. 

1.1 Alfven weak turbulence: the kinematics, the asymptotic clo- 
sure and some results 

Weak turbulence theory is a widely familiar approach to the plasma physics community, 
see e.g. Vedenov (1967, 1968), Sagdeev and Galeev (1969), Tsytovich (1970), Kuznetsov 
(1972, 1973), Zakharov (1974, 1984), Akhiezer et al. (1975), Mclvor (1977), Achterberg 
(1979) and Zakharov et al. (1992). This approach considers statistical states which can 
be viewed as large ensembles of weakly interacting waves and which can be described by a 
kinetic equation for the wave energy. Recall that the IK theory considers large ensembles of 
weakly interacting Alfven waves, but IK do not derive a kinetic equation and they restrict 
themselves to phenomenology based on the dimensional argument. Ng and Bhattacharjee 
(1996) developed a theory of weakly interacting Alfven wave packets which takes into ac- 
count anisotropy which leads to certain predictions for the turbulence spectra based on some 
additional phcnomenological assumptions and by-passing derivation of the weak turbulence 
kinetic equations. To date, there exists no rigorous theory of weak Alfven turbulence in 
incompressible MHD, and derivation of such a theory via a systematic asymptotic expansion 
in powers of small nonlinearity is the main goal of the present paper. It is interesting that the 
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kinetic equations were indeed derived (in some limits) for the Alfven waves for the cases when 
such effects, as finite Larmor radius (Mikhailovskii at al. 1989) or compressibihty (Kaburaki 
and Uchida 1971; Kuznetsov 1973), make these waves being dispersive. Perhaps the main 
reason why such a theory has not been developed for the non-dispersive Alfven waves in 
incompressible magneto- fluids was a general understanding within the "weak turbulence" 
community that a consistent asymptotic expansion is usually impossible for nondispersive 
waves. The physical reason for this is that all wavepackets propagate with the same group 
velocity even if their wavenumbers are different. Thus, no matter how weak the nonlinearity 
is, the energy exchanged between the wavepackets will be accumulated over a long time and 
it may not be considered small, as it would be required in the weak turbulence theory. As we 
will show in this paper, the Alfven waves represent a unique exception from this rule. This 
arises because the nonlinear interaction coefficient for co-propagating waves is null, whereas 
the counter-propagating wavepackets pass through each other in a finite time and exchange 
only small amounts of energy, which makes the weak turbulence approach applicable in this 
case. Because of this property, the theory of weak Alfven turbulence which is going to be 
developed in this paper posesses a novel and interesting mathematical structure which is 
quite different from the classical weak turbulence theory of dispersive waves. 

The starting point in our derivation is a kinematic description of the fields. We assume 
that the Elsasser fields z^(x, t) are random, homogeneous, zero mean fields in the three spatial 
coordinates x. This means that the n-point correlation functions between combinations of 
these variables estimated at xi, ...,Xn depend only on the relative geometry of the spatial 
configuration. We also assume that for large separation distances |xi — Xj| along any of the 
three-spatial directions, fluctuations are statistically independent. We will also discuss the 
case of strongly two dimensional flelds for which there is signiflcant correlation along the 
direction of the applied magnetic field. We choose to use cumulants rather than moments, 
to which the cumulants are related by a one-to-one map. The choice is made for two reasons. 
The first is that they are exactly those combinations of moments which are asymptotically 
zero for all large separations. Therefore they have well defined and, at least initially before 
long distance correlations can be built up by nonlinear couplings, smooth Fourier transforms. 
We will be particularly interested in the spectral densities 

1 r+oo 

= J^s (4(x)4(x + r))e-'''-dr (1) 

of the two point correlations. (Remember, 2;|(x) has zero mean so that the second order 
cumulants and moments are the same.) The second reason for the choice of cumulants as 
dependent variables is that, for joint Gaussian fields, all cumulants above second order are 
identically zero. Moreover, because of linear wave propagation, initial cumulants of order 
three and higher decay to zero in a time scale (6o^||)~^ where bo = Bq/^/XoPo is the Alfven 
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velocity (69 = |bo|) and k^^ a dominant parallel length scale in the initial field. This is a 
simple consequence of the Riemann-Lebesgue lemma ; all Fourier space cumulants become 
multiplied by fast nonvanishing oscillations because of linear wave properties and these os- 
cillations give rise to cancelations upon integration. Therefore, the statistics approaches a 
state of joint Gaussianity. The amount by which it differs, and the reason for a nontrivial 
relaxation of the dynamics, is determined by the long time cumulative response generated by 
nonlinear couplings of the waves. The special manner in which third and higher cumulants 
are regenerated by nonlinear processes leads to a natural asymptotic closure of the statistical 
initial value problem. 

Basically, because of the quadratic interactions, third order cumulants (equal to third 
order moments) are regenerated by fourth order cumulants and binary products of second 
order ones. But the only long time contributions arise from a subset of the second order 
products which lie on certain resonant manifold defined by zero divisors. It is exactly these 
terms which appear in the kinetic equations which describe the evolution of the power spec- 
tra of second order moments over time scales (e^ 6o^||)~^- Here, e is a measure of the strength 
of the nonlinear coupling. Likewise, higher order cumulants are nonlinearly regenerated 
principally by products of lower order cumulants. Some of these small divisor terms con- 
tribute to frequency renormalization and others contribute to further {e.g. four wave resonant 
interactions) corrections of the kinetic equations over longer times. 

What are the resonant manifolds for three wave interactions and, in particular, what 
are they for Alfven waves ? They are defined by the divisors of a system of weakly coupled 
wave-trains aje'^'^^'^'^''^^^^*\ with ujs{kj) the linear wave frequency, s its level of degeneracy, 
which undergo quadratic coupling. One finds that triads k, k, L which lie on the resonant 
manifold defined for some choice of s, s', s", by 

k = K -I- L , 

interact strongly (cumulatively) over long times (e^c^o)"^, ujq being a typical frequency. For 
Alfven waves, ujsO^) = sho ■ k = s6o^|| when s = ±1 (Alfven waves of a given wavevector 
can travel in one of two directions) and 60, the Alfven velocity, is the strength of the applied 
field. Given the dispersion relation, uj = sho ■ k, one might ask why there is any weak 
turbulence for Alfven waves at all because for, s = s' = s", (0) is satisfied for all triads. 
Furthermore, in that case, the fast oscillations multiplying the spectral cumulants of order 
iV + 1 in the evolution equation for the spectral cumulant of order disappear so that there 
is no cancelation (phase mixing) and therefore no natural asymptotic closure. However, the 
MHD wave equations have the property that the coupling coefficient for this interaction is 
identically zero and therefore the only interactions of importance occur between oppositely 
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(2) 



traveling waves where s' = — 



s, s" = s. In this case, becomes 



2sbo ■ K = 2sboK 



0. 



(3) 



The third wave in the triad interaction is a fluctuation with zero Alfven frequency. SG 
incorrectly conclude that the effective amplitude of this zero mode is zero and that therefore 
the resonant three wave interactions are null. 

Although some of the kinetic equations will involve principal value integral (PVI) with 
denominator sco'(k) + suj{k) — suj{h — k) = 2s6oK||, whose meaning we discuss later, the 
majority of the terms contain the Dirac delta functions of this quantity. The equation for 
the total energy density contains only the latter implying that energy exchange takes place 
by resonant interactions. Both the resonant delta functions and PVI arise from taking long 
time limits t — > oo, e^t finite, of integrals of the form 



Therefore, implicit in the derivation of the kinetic equations is the assumption that F{k\\, t^t) 



is the integral over of a product of a combination of gj|/(k_L — k\\) with Q^'^{k^, 0) = 
Sp m^p^m Q'pm~'^('^^5 0)- Three observations (01,2,3) and two questions (Ql,2) arise from this 
result. 

01- Unlike the cases for most systems of dispersive waves, the resonant manifolds for 
Alfven waves foliate wavevector space. For typical dispersion relations, a wavevector k, lying 
on the resonant manifold of the wavevector k, will itself have a different resonant manifold, 
and members of that resonant manifold will again have different resonant manifolds. Indeed 
the union of all such manifolds will fill k space so that energy exchange occur throughout 
all of k space. 

02- In contrast, for Alfven waves, the kinetic equations for the total energy density 
contains k\\ as a parameter which identifies which wavevector plane perpendicular to Bq we 
are on. Thus the resonant manifolds for all wavevectors of a given k\\ is the plane k\\ = k\\. 
The resonant manifolds foliate k-space. 

03- Further, conservation of total energy holds for each k\\ plane. There is energy ex- 
change between energy densities having the same k\\ value but not between those having 
different k\\ values. Therefore, relaxation towards a universal spectrum with constant trans- 
verse flux occurs in wavevector planes perpendicular to the applied magnetic field. The 





F{k\\,e^t) remains nearly constant for k\\ ~ e 



2. In 



(5) 
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dependence of the energy density on /cy is nonuniversal and is inherited from the initial 
distribution along k\\. 

Ql- If the kinetic equation describes the evolution of power spectra for values of k\\ 
outside of a band of order e^, ^ < 2, then how does one define the evolution of the quantities 
contained in Q^^{k±, 0) so as to close the system in k\\ ? 

Q2- Exactly what is Q~^{k,±,0) ? Could it be effectively zero as SG surmise? Could it 
be possibly singular with singular support located near k\\ = in which case the limit (^) is 
suspect ? 

To answer the crucially important question 2, we begin by considering the simpler exam- 
ple of a one dimensional, stationary random signal u{t) of zero mean. Its power spectrum 
is f{uj) the limit of the sequence fii^^) = ^ I-L{'^it)'^it + T))e~*'^'^ rfr which exists because 
the integrand decays to zero as r — > ±oo. Ergodicity and the stationarity of u{t) allows us 
to estimate the average -R(t) = {u{t)u{t + r)) by the biased estimator 

1 fL-\T\/2 

^Lir) = — u{t- T/2)u{t + r/2) dt 

IL J-L+\t\/2 

with mean E{Rl{t)} = (1 — t /2L)R(t). Taking L sufficiently large and assuming a suffi- 
ciently rapid decay so that we can take Rhi'^) = for |r| > 2L means that Rl{t) is simply 
the convolution of the signal with itself. Furthermore the Fourier transform Sl^uj) can then 
be evaluated as 

1 /■2L 1 rL 

Sl{uj) = — / i?T(r)e-^-^ dT = -—\ u{t)e-'^' dt\^ . 

271 J-2L A7rL J-L 

For sufficiently large L, the expected value of Sl^u) is S^u), the Fourier transform of R{t) 
although the variance of this estimate is large. Nevertheless Sl{uj), and in particular Sl{0), 
is generally non zero and measures the power in the low frequency modes. To make the 
connection with Fourier space, we can think of replacing the signal u{t) by the periodic 
extension of the truncated signal UL{t) = u{t), \t\ < L ; UL{t + 2L) for \t\ > L. The zero 
mode of the Fourier transform 0^(0) = ^ J^i^u{t) dt is a nonzero random variable and, 
while its expected value (for large L) is zero, the expected value of its square is certainly 
not zero. Indeed the expected value of Sl{0) = 2La|(0) has a finite nonzero value which, as 
L oo, is independent of L as 0^(0) has zero mean and a standard deviation proportional 
to (2L)~^/^. Likewise for Alfven waves, the power associated with the zero mode Q~^{k±, 0) 
is nonzero and furthermore, for the class of three dimensional fields in which correlations 
decay in all directions, Q~''^{k±, k\\) is smooth near k\\ = 0. Therefore, for these fields, we may 
consider Q~''^{ti±, 0) as a limit of Q~^{k±, k\\) as ^ and oo. Here fcj_o is some 

wavenumber near the energy containing part of the inertial range. Therefore, in this case, 
we solve first the nonlinear kinetic equation for limfcu^o e|(A;_L, k\\), namely for very oblique 
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Alfven waves, and having found the asymptotic time behavior of e^{k±,0), then return to 
solve the equation for e^{k±, k\\) for finite k\\. 

Assuming isotropy in the transverse k± plane, we find universal spectra c^/c"* for E^{k±) 
(/ E^{k±, 0) dk± — J e^{k±) dk±), corresponding to finite fluxes of energy from low to high 
transverse wavenumbers. Then e*(A;_L, /cy) = f^{k\\)Cnk±"''^ where /*(A;||) is not universal. 
These solutions each correspond to energy conservation. We find that convergence of all 
integrals is guaranteed for —3<n±<—l and that 

n+ + n_ = —4 (6) 

which means, that for no directional preference, n+ — — —2. These solutions have finite 
energy, i.e. J E dk± converges. If we interpret them as being set up by a constant flux of 
energy from a source at low k± to a sink at high k±, then, since they have finite capacity 
and can only absorb a finite amount of energy, they must be set up in finite time. When we 
searched numerically for the evolution of initial states to the final state, we found a remark- 
able result which we yet do not fully understand. Each E^{k±) behaves as a propagating 
front in the form E'(k±) = {to - ty^^Eo(k±(to - tf'^) and Eo{l) ~ l'^^ as I +oo. This 
means that for t < to, the E^{k±) spectrum had a tail for k± < {to — t)'^^'^ with stationary 
form k]^^^ joined to = through a front Eo(k±(to — i)^/^). The 7/3 spectrum is steeper 
than the +2 spectrum. Amazingly, as t approached very closely to to, disturbances in the 
high k± part of the k^ solution propagated back along the spectrum, rapidly turning it 
into the finite energy flux spectrum A;J^. Wc neither understand the origin nor the nature of 
this transition solution, nor do we understand the conservation law involved with the second 
equilibrium solution of the kinetic equations. Once the connection to infinity is made, how- 
ever the circuit between source and sink is closed and the finite flux energy spectrum takes 
over. 

To this point we have explained how MHD turbulent fields for which correlations decay 
in all directions relax to quasiunivcrsal spectra via the scattering of high frequency Alfven 
waves with very oblique, low frequency ones. But there is another class of fields that it is 
also important to consider. There are homogeneous, zero mean random fields which have 
the anisotropic property that correlations in the direction of applied magnetic field do not 
decay with increasing separation Bq ■ (xi — X2). For this case, we may think of decomposing 
the Elsasser fields as 

z]{^±,x\\) = ^j(xx) +i|(x^,.X||) (7) 

where the Zj{'x±,x\\) have the same properties of the fields considered heretofore but where 
the average of Zj{x±,x\\) over x\\ is nonzero. The total average of Zj is still zero when one 
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averages Zj{x.±) over x^. In this case, it is not hard to show that correlations 

{zj{x±, X||) 4'(xx + r±, X|| + r||)) 

divide into two parts 

{z'j{x±) z]',(x± + r±)) + (ij(xi, ij,'(xx + r^, + ry)) 
with Fourier transforms, 

q^}ik) = q]}{k^ + qg,{k^, , (8) 

when q is smooth in and q is smooth in both and k\\. The former is simply the 
transverse Fourier of the two point correlations of the x\\ averaged field. Likewise all higher 
order cumulants have delta function multipliers S{k\\) for each k dependence. For example 

q^}]:{k,k') = 6{k^6{k\^)q]}p{k^,k'^ 

is the Fourier transform of 

(^j(xi) zf,{x± + rj_) zfi{x± + r^)) . 

Such singular dependence of the Fourier space cumulants has a dramatic effect on the 
dynamics especially since the singularity is supported precisely on the resonant manifold. 
Indeed the hierarchy of cumulant equations for ^("^ simply loses the fast (Alfven) time 
dependence altogether and becomes fully nonlinear MHD turbulence in two dimensions with 
time t replaced by et. Let us imagine, then, that the initial fields are dominated by this 
two dimensional component and that the fields have relaxed on the time scale t ~ e"^ to 
their equilibrium solutions of finite energy flux for which E{k±) is the initial Kolmogorov 
finite energy flux spectrum kj^^^ for k± > ko, ko some input wavenumber and E(k±) ~ 
kj^^^ corresponding to the inverse fiux of the squared magnetic vector potential ; b = 
V X ^). ^(k), the spectral density of {A'^), behaves as kj^^^. These are predicted from 
phenomenological arguments and supported by numerical simulations. 

Let us then ask : how do Alfven waves (Bragg) scatter off this two dimensional turbulent 
field? To answer this question, one should of course redo all the analysis taking proper 
account of the factors in q^'"\ However, there is a simpler way. Let us imagine that the 
power spectra for the Zj fields are supported at finite k\\ and have much smaller integrated 
power over an interval < k\\ < (3 1 than do the two dimensional fields. Let us replace 
the S{k\\) multiplying gjj'(kx) by a function of finite width P and height Then the 

kinetic equation is linear and describes how the power spectra, and in particular e*(k). 
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of the Zj fields interact witfi the power spectra of the two dimensional field Zj. Namely, 
the Q~^{k±, 0) field in the kinetic equation is determined by the two dimensional field and 
taken as known. The time scale of the interaction is now /9e~^, because the strength of the 
interaction is increased by P'^ and, is faster than that of pure Alfven wave scattering. But 
the equilibrium of the kinetic equation will retain the property that n_s + Us — —4: where 
now n_s is the phenomenological exponent associated with two dimensional MHD turbulence 
and Ug the exponent of the Alfven waves. Note that when n_s = —5/3, is —7/3, which 
is the same exponent (perhaps accidentally) as for the temporary spectrum observed in the 
finite time transition to the fcj^ spectrum. 

We now proceed to a detailed presentation of our results. 

2 The derivation of the kinetic equations 

The purpose in this section is to obtain closed equations for the energy and helicity spectra 
of weak MHD turbulence, using the fact that, in the presence of a strong uniform magnetic 
field, only Alfven waves of opposite polarities propagating in opposite directions interact. 

2.1 The basic equations 

We will use the weak turbulence approach, the ideas of which are described in great detail 
in the book of Zakharov et al. (1992). There are several different ways to derive the weak 
turbulence kinetic equations. We follow here the technique that can be found in Benney and 
Newell (1969). We write the 3D incompressible MHD equations for the velocity v and the 
Alfven velocity b 

(9t + v- V)v = -VP, + b- Vb + i/V^v, (9) 

(^t + v- V)b = b- Vv + r^V^b, (10) 

where is the total pressure, u the viscosity, 1] the magnetic diffusivity and V ■ v = 0, 
V • b = 0. In the absence of dissipation, these equations have three quadratic invariants in 
dimension three, namely the total energy = + the cross-correlation E'-^ — (v-b) 
and the mag netic helicity = ■ b) (Woltjer 1958). 

The Elsasser variables = v + sb with s = ±1 give these equations a more symmetrized 
form, namely : 

{dt + z-'-V)z' = -VP. , (11) 

where we have dropped the dissipative terms which pose no particular closure problems. The 
first two invariants are then simply written as 2E^ — (|z*p). 
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We now assume that there is a strong uniform magnetic induction field Bq along the unit 
vector ey and non dimensionalize the equations with the corresponding magnetic induction 
Bo, where the z'^ fields have an amplitude proportional to e (e <^ 1) assumed small compared 
to 60 • Linearizing the equations leads to 

{dt - s6o9||)^j = -e9.„^-^^| - , (12) 

where d\\ is the derivative along ey. The frequency of the modes at a wavevector k is 
co'(k) = cjfc = bo ■ k = 60^11 • We Fourier transform the wave fields using the interaction 
representation, 

;2|(x, t) = j Aj(k, t) e**^ " rfk = y aj(k, t) e^(k-+--fet) d\<. , (13) 

where a*(k, t) varies slowly in time because of the weak nonlinearities ; hence 

dta]{}^, t) = -zek^P.n J a-\K) <(L) g^C— (14) 

with c/kl = dndL and 5k,KL = ^(k — k — L) ; finally, Pjn{k) = 6jn — kjknk~^ is the usual 
projection operator keeping the A^(k) fields transverse to k because of incompressibility. 
The exponentially oscillating term in (|14|) is essential : its exponent should not vanish when 
(k — K — L) = 0, i.e. the waves should be dispersive for the closure procedure to work. 
In that sense, incompressible MHD can be coined "pseudo" -dispersive because although 
LJk ~ k, the fact that waves of one s-polarity interact only with the opposite polarity has 
the consequence that the oscillating factor is non-zero except at resonance ; indeed with 
ujic = bokw, one immediately sees that —suk — su^ + su^ = s{—k\\ — /ty + L\\) = —2sk,\\ using 
the convolution constraint between the three waves in interaction. In fact, Alfven waves 
may have a particularly weak form of interactions since such interactions take place only 
when two waves propagating in opposite directions along the lines of the uniform magnetic 
field meet. As will be seen later (see §|]), this has the consequence that the transfer in the 
direction parallel to Bq is zero, rendering the dynamics two-dimensional, as is well known 
(see e.g. Montgomery and Turner 1981; Shebalin et al. 1983). Technically, we note that there 
are two types of waves that propagate in opposite directions, so that the classical criterion 
(Zakharov et al. 1992) for resonance to occur, viz. u" > does not apply here. 

2.2 Toroidal and poloidal fields 

The divergence-free condition implies that only two scalar fields are needed to describe the 
dynamics ; following classical works in anisotropic turbulence, they are taken as (Craya 1958; 
Herring 1974; Riley et al. 1981) 

= + = V X (^^e||) + V X (V X (^^e,,)) , (15) 
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which in Fourier space gives 

A'j{k) = ik X e|| ^"^(k) - k x (k x ey) 0^(k) . (16) 

We elaborate somewhat on the significance of the ip'^ and (p^ fields since they are the basic 
fields with which we shall deal. Note that are two-dimensional fields with no parallel 
component and thus with only a vertical vorticity component (vertical means parallel to 
Bo), whereas the fields have zero vertical vorticity; such a decomposition is used as well 
for stratified fiows (see Lesieur 1990 and references therein). Indeed, rewriting the double 
cross product in (|16]) leads to : 

A'(k) = ikx ey ^^(k) -kk\\(f)'{k) + ey k'^(f)'{k) (17) 

or using k = k_L + A;||e|| : 

A"(k) = ikx e|| ^^^(k) - k^k^ 0*(k) + ey kl4)'{k) . (18) 

The above equations indicate the relationships between the two orthogonal systems (with 
p = k X e|| and q = k x p) made of the triads (k, p, q), (ey, p, k_|_) and the system (ey, p, k). 
In terms of the decomposition used in Waleffe (1992) with 

h± = pxk±zp (19) 

and writing z* = y4^h++A!_h_, it can be shown easily that = A^_^ — At and 0* = A^ + Ai. 
In these latter variables, the s-energies E''' are proportional to + and the s- 

helicities (z'' ■ V x z'^) are proportional to — Note that E'"^ is not a scalar: 

when going from a right-handed to a left-handed frame of reference, E"^ changes into E~^. 

2.3 Moments and cumulants 

We now seek a closure for the energy tensor qjj'ik) defined as 

(a^^(k)a^<(k')) = glf:(k')<5(k + k') (20) 

in terms of second order moments of the two scalar fields ip^{k) and 0*(k). Simple manipu- 
lations lead, with the restriction s = s' (it can be shown that correlations with s' = —s have 
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no long time influence and therefore are, for convenience of exposition, omitted), to : 

qttik') = A;2^^(k) - kik2kiir{k) + A;jA;2$^(k) , 

9i2(k') + (k') = -2A;iA;2*'(k) + A;||(A;2 - k:l)P{k) + 2A;iA;2A;J$^(k) , 

qllik') + qllik') = k2klr{k)-2k,k\\kl<^^{k), ^^^^ 

g|f (k') + gp(k') = -k^klP{k) - 2k2k^^kl^^ (k) , 

gfll(k') = ^l$^(k), 

-b2l(k') - gp(k')] = i-Mli^') - 5if (k')] = i-Mtm - 5l?(kO] 



fcl 



-iA;ii?*(k) , 



where the following correlators involving the toroidal and poloidal fields have been intro- 
duced : 

{^'{k)^'{k')) = S{k + k')^'(k') , 

(0-(k)0«(k')) = S{k + k')$^(k') , 

(z/J«(k)0-(k')) = S{k + k')n^(-k) , 

. . (22) 
(0-(k)^-(k')) =5(k + k')ff(k), 

R'{k) = n*(-k) + n*(k) , 

7^(k) = i[n^(-k) - n^(k)] , 

and where k'j_ = kj + k^ , k'^ = kj_ + fcjj. Note that E^i?* is the only pseudo-scalar, linked to 
the lack of symmetry of the equations under plane reversal, i.e. to a non-zero helicity. 
The density energy spectrum writes 

e^(k) = Ej gjj(k) = kl(*^(k) + k^^'ik)) . (23) 

Note that it can be shown easily that the kinetic and magnetic energies and |(&^) 

are equal in the context of the weak turbulence approximation. Similarly, expressing the 
magnetic induction as a combination of z"^ and thus of ■^^ and 0^, the following symmetrized 
cross-correlator of magnetic hehcity (where the Alfven velocity is used for convenience) and 
its Fourier transform are found to be 

^(A(k) 6,(k')) + liMk') 6,(k)) = ^kiE,i?^(k) 6{k + k') , (24) 

where the correlations between the + and — variables are ignored because they are expo- 
nentially damped in the approximation of weak turbulence. Similarly to the case of energy. 
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there is equivalence between the kinetic and magnetic hehcal variables in that approxima- 
tion, hence the kinetic helicity defined as (u ■ a;) writes simply in terms of its spectral density 
H^{\.) : 

H^{\i) = k^H^iy) = h^lil^sR'iy) . (25) 

In summary, the eight fundamental spectral density variables for which we seek a weak 
turbulence closure are the energy e^(k) of the three components of the fields, the energy 
density along the direction of the uniform magnetic field $*(k), the correlators related to the 
off-diagonal terms of the spectral energy density tensor /*(k) and finally the only helicity- 
related pseudo-scalar correlators, namely i?*(k). 

The main procedure that leads to a closure of weak turbulence for incompressible MHD 
is outlined in the Appendix. It leads to the equations (|110|) giving the temporal evolution of 



the components of the spectral tensor q'j|''(k) just defined. The last technical step consists 
in transforming equations ( |110| ) of the Appendix in terms of the eight correlators we defined 
above. This leads us to the final set of equations, constituting the kinetic equations for weak 
MHD turbulence. 

2.4 The kinetic equations 

In the general case the kinetic equations for weak MHD turbulence are 



(26) 



Tie 



fci-^W^(k) + UiL^-^J$W 



klY^' 



2 l2 



}\i I k 



klX^ 



dt [k\k'^^\\i) 
'^^(L) $^(k)\ 



(27) 



klk^ 



+Mr {qz + klLl) r (L) + ( MZ") r(k) 



2L2 



2 /^"(L) <l>^(k)\ 



^^^^(k)P / 



X 
2kiiL2 



14 



klR'{k) 



(28) 



ire 



'Z + kl 



A;2 / {Z + kh^\ 



Qr(«) 



+ (i^k\\Z-L\^kiy -ex^^ r(k) 



with 



and 



a 



A;lA;'7^(k) 



Tie 



LlZ+^iZ'-X')\nL) + 



r(k) 



2A;iiX 



2L2 ^ 2L2 ^ 
+ (^^^) (*1k) + ^2$^(k)) + ^(Z*^(L) - [kp + klLl) $^(L)) 



(^k,icL = 5(L + k;- k), 



(29) 



(30) 



= X2*-^(«) + X(A;||Ki - K||r)/-*(«) + {kiiY - A;||/«i) V(«) ■ 

Note that does not involve the spectral densities R^(k.), because of symmetry prop 
erties of the equations. The geometrical coefficients appearing in the kinetic equations are 



X — (kj_ A k_l)2 = A;_L«;_L sin 6* , 

Y — Icj^ - K± — k±K,± cos 9 , 

Z — 'k±-L± — kj_ — k±K± cos 9 — k^ — Y , 

W — K± ■ L± = k'j_ — L\ — k±K± cos 9 — Z — L\, 



(31) 
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where 6 is the angle between and k^, and with 



dK± = n±dK±d9 



cos 9 



sin 9 

Kj_ -|- k'j_ -Zyj_ 



dK\dL\ 



(32) 
(33) 



In (^Tf), (pSD and (E^) P / means the Cauchy Principal value of the integral in question. 



3 General properties of the kinetic equations 

3.1 Dynamical decoupling in the direction parallel to Bo 

The integral on the right-hand side of the kinetic equation ( p6D contains a delta function of 
the form the integration variable corresponding to the parallel component of one of the 

wavenumbers in the interacting triad. This delta function arises because of the three-wave 
frequency resonance condition. Thus, in any resonantly interacting wave triad (k, k, L), 
there is always one wave that corresponds to a purely 2D motion - having no dependence 
on the direction parallel to the uniform magnetic field - whereas the other two waves have 
equal parallel components of their corresponding wavenumbers, viz. L\\ = k\\. This means 
that the parallel components of the wavenumber enter in the kinetic equation of the total 
energy e*(k) as an external parameter and that the dynamics is decoupled at each level of k\\. 
In other words, there is no transfer associated with the three-wave resonant interaction along 
the /c||-direction in k-space for the total energy. This result, using the exact kinetic equations 
developed here, corroborates what has already been found in Montgomery and Turner (1981) 
using a phenomenological analysis of the basic MHD equations, in Ng and Bhattacharjee 
(1996, 1997) in the framework of a model of weak MHD turbulence using individual wave 
packets, and in Kinney and McWilliams (1998) with a Reduced MHD approach (RMHD). 

As for the kinetic equation (^), the other kinetic equations ( pTj) to (p9D have integrals 
containing delta functions of the form S{k,\\). But, in addition, they have PVIs which can, a 
priori, contribute to a transfer in the parallel direction. The eventual contributions of these 
PVIs are discussed in § p.4| . 



3.2 Detailed energy conservation 

Detailed conservation of energy for each interacting triad of waves is a usual property in weak 
turbulence theory. This property is closely related with the frequency resonance condition 
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because uj can be interpreted as the energy of one wave "quantum". For Alfven waves, 
the detailed energy conservation property is even stronger because one of the waves in any 
resonant triad belongs to the 2D state with frequency equal to zero, 

ojii oc K|| = 0. 

Thus, for every triad of Alfven waves k, L and k (such that n\\ = 0) the energy is conserved 
within two co-propagating waves having wavevectors k and L. Mathematically, this corre- 
sponds to the symmetry of the integrand in the equation for with respect to changing 
k ^ L (and correspondingly k = k — L ^ —k). 

As we have said, energy is conserved fcy plane by fcy plane so that, for each fcy, it can be 
shown from (p^) 

^y'e^(k^,fc||)rfk^ = 0. (34) 
3.3 The magnetic and pseudo magnetic helicities 

Since Woltjer (1958) we know that the magnetic helicity is an invariant of the MHD equa- 
tions. However Stribling et al. (1994) showed that in presence of a mean magnetic field 
Bo the part of the magnetic helicity associated with fluctuations is not conserved separately 
(whereas the total magnetic helicity, which takes into account a term proportional to Bq, is, 
of course, an invariant). It is then interesting to know if, in the context of weak turbulence, 
the integral of the spectral density of fluctuations of the magnetic helicity is conserved, i.e. 

J H^{k) dk = constant , (35) 

with 

H^\k) = h^klE.R'ik) . (36) 
To investigate this point we define in the physical space the total magnetic helicity as 

= {At ■ bT) , (37) 

where bx = V x Ar and bx = bo + b. The magnetic induction equation 

dthT = V X (v X bx) (38) 

implies that (Stribling et al. 1994) 

dtH^' = dtH^' + 2ho-dt{A), (39) 

where is the magnetic helicity associated with fluctuations {H^ = ■ b)) and b = 
V X Direct numerical simulations (Stribling et al. 1994) show that the second term in the 
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RHS of ( p9D has a non-zero contribution to the total magnetic hehcity, but in the context 
of weak turbulence the situation is different. Indeed, the magnetic induction equation leads 
also to the relation 

9,(^) = ^(z-xz+). (40) 

Therefore the temporal evolution of the magnetic potential of fluctuations is proportional to 
the cross product between z- fields of opposite polarities. As we have already pointed out, 
in the framework of weak turbulence this kind of correlation has no long time influence and 
thus the magnetic helicity associated with fluctuations appears to be an invariant of the weak 
turbulence equations. We leave for the future the investigation of this point : in particular it 
will be helpful to make numerical computations to show, at least at this level, the invariance 
of the magnetic helicity. 

The correlators R^(k) and I^(k) have been defined in the previous section as the real 
part and the imaginary part of 11* (k), the cross-correlator of the toroidal field '0*(k) and of 
the poloidal field 0*(k). Then k5^Es-R''(k) appears as the spectral density of the magnetic 
helicity. On the other hand /*(k), which we will call the anisotropy correlator (or pseudo 
magnetic helicity), is neither a conserved quantity nor a positive definite quantity. Although 

and J* evolve according to their own kinetic equations ( pOj) and (pSj), the range of values 
they can take on is bounded by \Ef* and with the bounds being a simple consequence of 
the definition of these quantities. Two realizability conditions (see also Cambon and Jacquin 
1989; Cambon et al. 1997) between the four correlators \E'*, P and can be obtained 
from 

{\i!'{k)±k(f)'{k)\^) >0, (41) 

and 

(if(k)n(i<^"^(k)n>i(^^(k)0^(-k))p. (42) 

These conditions are found to be respectively 

**(k) + k^^'ik) > \kR'{k)\ , (43) 

and 

4^"(k)$*(k) > R'\k) + r\k) . (44) 

Note that the combination 

Z = {l/2)kl[e^k) - ^(k) - t\k\I{k)] (45) 

is named polarization anisotropy in Cambon and Jacquin (1989). The consequences of the 
realizability conditions is explained below. 
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3.4 Purely 2D modes and two-dimensionalisation of 3D spectra 

The first consequence of tlie fact that there is no transfer of the total energy in the /cy 
direction in k-space is an asymptotic two-dimensionahsation of the energy spectrum e*(k). 
Namely, the 3D initial spectrum spreads over the transverse wavenumbers, k_|_, but remains 
of the same size in the fcy direction, and the support of the spectrum becomes very fiat 
(pancake-like) for large time. The two-dimensionalisation of weak MHD turbulence has 
been observed in laboratory experiments (Robinson and Rusbridge 1971), in the solar wind 
data (Bavassano et al. 1982; Matthaeus et al. 1990; Horbury et al. 1995; Bieber et al. 
1996), and in many direct numerical simulations of the three-dimensional MHD equations 
(Oughton et al. 1994) or of the RMHD equations (Kinney and McWilliams 1998). 

From a mathematical point of view, the two-dimensionalisation of the total energy means 
that, for large time, the energy spectrum e''(k) is supported on a volume of wavenumbers 
such that for most of them k± 3> fcy. This implies that \E''^(k) and $'^(k) are also supported 
on the same anisotropic region of wavenumbers because both of them are non-negative. This, 
in turn, implies that both and P will also be non-zero only for the same region in the 
k-space as e'^(k), \E''^(k) and $^(k), as it follows from the bound (^) and d^). This fact 
allows one to expand the integrands in the kinetic equations in powers of small k\\/k±. At 
the leading order in k\\/k±, one obtains 



Li - 




V2\ 

— r(k) 



(46) 



ne 



'0 



dt[ki<!>%k)] 



(47) 



dtlklR^m 




(4J 




(49) 
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Note that the principal value terms drop out of the kinetic equations at leading order. This 
property means that there is no transfer of any of the eight correlators in the direction in 
k-space. 

One can see from the above that the equations for the toroidal and poloidal energies 
decouple for large time. These equations describe the shear-Alfven and pseudo-Alfven waves 
respectively. An energy exchange between \E'**(k) and $'*(k) is however possible in an initial 
phase, i.e. before the two-dimensionalisation of the spectra. A preliminary investigation 
shows that this exchange is actually essentially generated by the magnetic helicity through 
the principal value terms : the magnetic helicity plays the role of a catalyst which transfers 
toroidal energy into poloidal energy. On the other hand, in the large time limit, the magnetic 
helicity T.gR'^ and the pseudo magnetic helicity P are also described by equations which are 
decoupled from each other and from the toroidal and poloidal energies. It is interesting 
that the kinetic equation for the shear-Alfven waves {i.e. for \E'*(k)) can be obtained also 
from the RMHD equations which have been derived under the same conditions of quasi 
two-dimensionality (see e.g. Strauss 1976). 

An important consequence of the dynamical decoupling at different /cy's within the kinetic 
equation formalism is that the set of purely 2D modes (corresponding to fcy = 0) evolve 
independently of the 3D part of the spectrum (with k\\ ^ 0) and can be studied separately. 
One can interpret this fact as a neutral stability of the purely 2D state with respect to 3D 
perturbations. As we mentioned in the Introduction, the kinetic equations themselves are 
applicable to a description of fcy = modes only if the correlations of the dynamical fields 
decay in all directions, so that their spectra are sufficiently smooth for all wavenumbers 
including the ones with k\\ =0. To be precise, the characteristic k\\ over which the spectra 
can experience significant changes must be greater than e^. Study of such 2D limits of 3D 
spectrum will be presented in the next section. It is possible, however, that in some physical 
situations the correlations decay slowly along the magnetic field due to a (hypothetical) 
energy condensation at the = modes. In this case, the modes with k\\ = should be 
treated as a separate component, a condensate, which modifies the dynamics of the 3D modes 
in a manner somewhat similar to the superfiuid condensate, as described by Bogoliubov 
(Landau and Lifshitz 1968). We leave this problem for future study. 

3.5 Asymptotic solution of the 3D kinetic equations 

The parallel wavenumber k\\ enters equations (^6D-(^9|) only as an external parameter. In 
other words, the wavenumber space is foliated into the dynamically decoupled planes k\\ = 0. 
Thus, the large-time asymptotic solution can be found in the following form. 
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vl>^(k^,fc||) = f,m)^^{k^,0), (50) 

$^(k^,fc||) = /2(fc||)$^(k^,0), (51) 

i?^(k^,fc||) = /3(A;|,)i?^(kx,0), (52) 

I'{^±,h) = /4(fc|,)r(k^,0), (53) 

where /«, (i = 1, 2, 3, 4) are some arbitrary functions of fcy satisfying the conditions fi{0) = 1 
(and such that the bounds (|i3| ) and ( ^^ are satisfied). Substituting these formulae into 
(^)- (|49|) , one can readily see that the functions fi drop out of the problem, and the solution 
of the 3D equations is reduced to solving a 2D problem for \E'*(k^, 0), $*'(k^, 0), -R*(k^, 0) 
and /'*(k^,0), which will be described in the next section. 



4 Two— dimensional problem 



Let us consider Alfven wave turbulence which is axially symmetric with respect to the ex- 
ternal magnetic field. Then P{k±,0) = because of the condition /*(— k) = — /'*(k). In 
the following, we will consider only solutions with = 0. (One can easily see that will 
remain zero if it is zero initially.) The remaining equations to be solved are 



dEl{k^,0) 



dt 



2 7 

[{eL-ekfstne Sr(«x,0) 



k^El{L^,0) - L^El{k^,0) 



dn\dL\ 



dEf,{k^,0) 



Txe 



'0 



sine — ET'\ 



dt 



k±EI{L^,0) - L^EI{k^,0) 



dn I dL 



where and are the unit vectors along k^ and respectively and 

El{k^,0) = kl^'{k^,0) , 
E^Jk^,0) = kl<l>%k^,0) , 



(54) 



(55) 



(56) 
(57) 



are the horizontal and the vertical components of the energy density. Thus, we reduced the 
original 3D problem to finding solution for the purely 2D state. It may seem unusual that 
strongly turbulent 2D vortices (no waves for k\\ = !) are described by the kinetic equations 
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obtained for weakly turbulent waves. Implicitly, this fact relies on continuity of the 3D 
spectra near k\\ = 0, so that one could take the limit k\\ 0. In real physical situations such 
continuity results from the fact that the external magnetic field is not perfectly unidirectional 
and, therefore, there is a natural smoothing of the spectrum over a small range of angles. 

The equation (|5^) corresponds to the evolution of the shear-Alfven waves for which the 
energy fluctuations are transverse to Bq whereas equation (|^) describes the pseudo-Alfven 
waves for which the fluctuations are along Bq. Both waves propagate along Bq at the same 



Alfven speed. Equation (|5J) describes the interaction between two shear-Alfven waves, 
E^, propagating in opposite directions. On the other hand, the evolution of the pseudo- 
Alfven waves depend on their interactions with the shear-Alfven waves. The detailed energy 
conservation of the equation (^ID implies that there is no exchange of energy between the 
two different kinds of waves. The physical picture in this case is that the shear-Alfven waves 
interact only among themselves and evolve independently of the pseudo-Alfven waves. The 
pseudo-Alfven waves scatter from the shear-Alfven waves without amplification or damping 
and they do not interact with each other. 

Using a standard two-point closure of turbulence (see e.g. Lesieur 1990) in which the 
characteristic time of transfer of energy is assumed known and written a priori, namely the 
EDQNM closure, Goldreich and Sridhar (1995) derived a variant of the kinetic equation ( |MD 
but for strong anisotropic MHD turbulence. In their analysis, the ensuing energy spectrum, 
which depends (as it is well known) on the phenomenological evaluation of the characteristic 
transfer time, thus differs from our result where the dynamics is self-consistent, closure being 
obtained through the assumption of weak turbulence. 

It can be easily verified that the geometrical coefficient appearing in the closure equation 
in Goldreich and Sridhar (1995) is identical to the one we find for the E^{h±, k\\) spectrum 
in the two-dimensional case. However, the two formulations, beyond the above discussion on 
characteristic time scales, differ in a number of ways : (i) We choose to let the flow variables to 
be non mirror-symmetric, whereas helicity is not taken into account in Goldreich and Sridhar 
(1995) where they have implicitly assumed = 0; (ii) However, because of the anisotropy 
introduced by the presence of a uniform magnetic field, one must take into account the 
coupled dynamics of the energy of the shear Alfven waves, the pseudo-Alfven wave and the 
pseudo magnetic helicity J* ; indeed, even if initially /'^ = 0, it is produced by wave coupling 
and is part of the dynamics, (iii) In three dimensions, all geometrical coefficients that depend 
on fc^ = fc^ + A;| have a fcy-dependence which is a function of initial conditions and again is 
part of the dynamics. 



4.1 Kolmogorov spectra 
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4.1.1 The Zakharov transformation 



The symmetry of the previous equations allows us to perform a conformal transformation, 
called the Zakharov transformation (also used in modeling of strong turbulence, see Kraich- 
nan 1967), in order to find the exact stationary solutions of the kinetic equations as power 
laws (Zakharov et al. 1992). This operation (see Figure |lD consists of writing the kinetic 
equations in dimensionless variables ui = K±/k± and uj2 = L±/k±, setting by k"^, and 
then rearranging the collision integral by the transformation 

oo[ = -, (58) 

UJ2 

and 

u', = — . (59) 

The new form of the collision integral, resulting from the summation of the integrand in its 
primary form and after the Zakharov transformation, is 

^'^^ J [ 2C02 ) V'^^^i ' 

{uj^'~^ - 1)(1 - 002'''^'"'''^) duji duJ2 . 

The collision integral can be null for specific values of n±. The exact solutions, called the 
Kolmogorov spectra, correspond to these values which satisfy 

n+ + n_ = -4 . (60) 

It is important to understand that the Zakharov transformation is not an identity transfor- 
mation, and it can lead to spurious solutions. The necessary and sufficient condition for a 
spectrum obtained by the Zakharov transformation to be a solution of the kinetic equation is 
that the right hand side integral in ( ^41) {i.e. before the Zakharov transformation) equation 
converges. This condition is called the locality of the spectrum and leads to the following 
restriction on the spectral indices in our case : 

-3<n±<-l. (61) 

A detailed study of the Kolmogorov spectrum locality will be given in section ||. 

In the particular case of a zero cross-correlation one has = E2 = E±_ ~ /c" with only 
one solution 

n = -2 . 

Note that the thermodynamic equilibrium, corresponding to the equipartition state for which 
the ffux of energy is zero instead of being finite as in the above spectral forms, corresponds 
to the choices n_|_ = n_ = 1 for both the perpendicular and the parallel components of the 
energy. 
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4.1.2 The Kolmogorov constants CK{ns) and C^(ns) 

The final expression of ttie Kolmogorov-like spectra found above as a function of tiie Kol- 
mogorov constant (generalised to MHD) CK{ns) and of the flux of energy Pl{k±) can be 
obtained in the following way. For a better understanding, the demonstration will be done 
in the simplified case of a zero cross-correlation. The generalization to the correlated case 
{E~^ 7^ E^) is straightforward. Using the definition of the flux, 

dtE4k±,0) ^ -dk^P^ki,) , (62) 

one can write the flux of energy as a function of the collision integral (with the new form 
of the integrand) depending on n. Then the limit n — > —2 is taken in order to have a 
constant flux P± with no more dependence in k±, as it is expected for a stationary spectrum 
in the inertial range. Here we have considered an inflnite inertial range to use the Zakharov 
transformation. Whereas the collision integral tends to zero when n —2, the limit with 
which we are concerned is not zero because of the presence of a denominator proportional 
to 2n + 4, and which is a signature of the dimension in wavenumber of the flux. Finally the 
"L'Hospital's rule" gives the value of P± from which it is possible to write the Kolmogorov 
spectrum of the shear- Alfven waves 

E^k^, 0) = P]/' CKi-2) kz' , (63) 

with the Kolmogorov constant 



and with the following form for the integral Ji (n) 



^' Jx=i Jy=-i (x-y)"+6 x + y)2-n ^ ^ 



2^-" -{x + yy-""] In dxdy . 



As expected, the calculation gives a negative value for the integral Ji(n) and for the particular 
value n — —2, we obtain Ck{—'2) — 0.585. Note that the integral Ji(n) converges only for 
-3 < n< -1. 

The generahzation to the case of non-zero cross-correlation gives the relations 

Et{k±, 0) EZ{k^, 0) = P+ CUris) kl" = C|(-n, - 4) kz" 



^/pfPZCKiris) CKi-ris - 4) k^^ , (66) 
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where the second formulation is useful to show the symmetry with respect to s. The com- 
putation of the Kolmogorov constant Ck as a function of —Ug is given in Figure ^ An 
asymmetric form is observed which means that the ratio P^/PZ is not constant, as we can 
see in Figure ^ where we plot this ratio as a function of — n^. We see that for any ratio 
P^ / P± there corresponds a unique value of n^, between the singular ratios P±/P± = +00 
for Hg = —3 and P±/P± = for = — 1. Thus, a larger flux of energy P^ corresponds to 
a steeper slope of the energy spectra E^{k±, 0) in agreement with the physical image that a 
larger flux of energy implies a faster energy cascade. 

In the zero cross-correlation similar demonstration for the pseudo-Alfven waves 

E\\{k±,0) leads to the relation 

E||(A;^,0) = P||Pr'/' C'j,{-2) kf , (67) 

with the general form of the Kolmogorov constant 



CM 

where the integral J2{n) is 



-2boJi{n) 



\ 7re2J2(n)J2(-n-4) 



(68) 



,0 r+°° /■! a/(x2 - 1)(1 - y^) 
Jx=i Jy=~i{x-y)''+^ {x + yy-"" ^ ' 

>^~" - (x + y)^'"] In (^^) dxdy . 

Note that the integral Jiin) converges only for —3 < n < — 1. The presence of the flux Pj_ 
in the Kolmogorov spectrum is linked to the presence of in the kinetic equation of E\\. 
A numerical evaluation gives C^(— 2) ~ 0.0675 whereas the generalization for the non-zero 
cross-correlation is 

E+(A;^,0)E|7(A;^,0) = A;!^ = 5lA _ 4) 

C'M) CA-ris - ^) kj_\ (70) 



PtP^ 



where the last formulation shows the symmetry with respect to s. The power laws of the 
spectra E^^ have the same indices than those of E^ and the Kolmogorov constant C is in 
fact related to C by the relation 

CAns) ^ Cx(-n.-4) 
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Therefore the choice of the ratio / P± determines not only Cxins) but also Cj^(ns), 
allowing for free choices of the dissipative rates of energy P^. 

The result of the numerical evaluation of C^(ns) is shown in Figure ^. An asymmetrical 
form is also visible; notice also that the values of C^(ns) {i.e. the constant in front of the 
parallel energy spectra) are smaller by an order of magnitude than those of Cxing) for the 
perpendicular spectra. 

4.2 Temporal evolution of the kinetic equations 
4.2.1 Numerical method 



Equations (0) and (p5D can be integrated numerically with a standard method, as for 
example presented in Leith and Kraichnan (1972). Since the energy spectrum varies smoothly 
with k, it is convenient to use a logarithmic subdivision of the k axis 

k = Sk T'^ , (72) 

where z is a non-negative integer ; Sk is the minimum wave number in the computation and 
F is the number of wave numbers per octave. F defines the refinement of the "grid" , and in 
particular it is easily seen that a given value of F introduces a cut-off in the degree of non- 
locality of the nonlinear interactions included in the numerical computation of the kinetic 
equations. But since the solutions are local, a moderate value of F can be used (namely, 
we take F = 4). Tests have nevertheless been performed with F = 8 and we show that no 
significant changes occur in the results to be described below. 

This technique allows us to reach Reynolds numbers much greater than in direct numerical 
simulations. In order to regularize the equations at large k, we have introduced dissipative 
terms which were omitted in the derivation of the kinetic equations. We take the magnetic 
Prandtl number [i'/t]) to be unity. For example, with 6k = 2~^, F = 8, imax = 225; this 
corresponds to a ratio of scales 2^^/2~^. Taking a wave energy Uq and an integral scale Lq 
both of order one initially, and a kinematic viscosity of u = 3.3 x 10~^, the Reynolds number 
of such a computation is Tie = Uq Lq/u ~ 10^. All numerical simulations to be reported 
here have been computed on an Alpha Server 8200 located at the Observatoire de la Cote 
d'Azur (SIVAM). 



4.2.2 Shear-AIfven waves 

In this paper, we only consider decaying turbulence. As a first numerical simulation we 
have integrated the equation ( |MD in the zero cross-correlation case {E^ = E~) and without 
forcing. Figure ^ (top) shows the temporal evolution of the total energy E±(t) with by 
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definition 

E^it) = / E^ik^,0)dk± , (73) 

where kmin and kmax have the values given in the previous section. The total energy is 
conserved up to a time to ~ 1.55 after which it decreases because of the dissipative effects 
linked to mode coupling, whereas the enstrophy / k'^E±{k) ci^k increases sharply (bottom of 
Figure The energy spectra at different times are displayed in Figure |^. As we approach 
the time to, the spectra spread out to reach the smallest scales {i.e. the largest wavenumbers) . 
For t > to, constant energy flux spectrum kj^ is obtained (indicated by the straight line). For 
times t significantly greater than to, we have a self-similar energy decay, in what constitutes 
the turbulent regime. 



4.2.3 Shear-AIfven versus pseudo-Alfven waves 

In a second numerical computation we have studied the system ([5^ ) (l55|) with an initial 
normalised cross-correlation of 80%. The following parameters have been used: 6k = 2~^, 
F = 4, imax = 105 and u = 6.4 10~^. Figure |^ (top) shows the temporal evolution of energies 
for the four different waves {E^ and Ef). The same behavior as that of Figure ^ (top) 
is observed, with a conservation of energy up to the time t^ ^ 5, and a decay afterwards; 
this decay is nevertheless substantially weaker than when the correlation is zero, since in the 
presence of a significant amount of correlation between the velocity and the magnetic field, it 
is easily seen from the primitive MHD equations that the nonlinearities are strongly reduced. 
On the other hand the temporal evolution of enstrophies (bottom) displays that the maxima 
for these four types of waves are reached at different times : the pseudo-Alfven waves are 
the fastest to reach their maxima at t ~ 5.5 vs. t ~ 7.5 for the shear-Alfven waves. Figure 
^ corresponds to the temporal evolution of another conserved quantity, the cross correlation 
Pa; defined as 

" ^ irf ■ 

where x symbolizes either ± or ||. As expected, is constant during an initial period (till 
t = t'o) and then tends asymptotically to one, but in a faster way for the pseudo-Alfven 
waves. This growth of correlation is well documented in the isotropic case (Matthaeus and 
Montgomery 1980) and is seen to hold as well here in the weak turbulence regime. Figures ^ 
and give the compensated spectra E^E2k\ and E^^^E^k^ respectively at different times. 
In both cases, from t = Q onward, a plateau is observed over almost four decades and remains 
flat for long times; this illustrates nicely the theoretical predictions (|6^) and (|70|) . 
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5 Front propagation 



The numerical study of the transition between the initial state and the final state, where the 
/cj^-spectrum is reached, shows two remarkable properties illustrated by Figure 11 and 12 



We show in Figure |TT] (top), in lin-log coordinates, the progression with time of the front 
of energy propagating to small scales ; more precisely, we give the wavenumber at time t 
with an energy of, respectively, 10~^^ (dash-dot line) and 10~^^ (solid line). Note that all 
curves display an abrupt change at to — 1-55, after which the growth is considerably slowed 
down. Using this data. Figure [Tl| (bottom) gives log(A;_L) as a function of log(1.55 — t), the 
lines having the same meaning as in Figure [ll| (top); the large dash represents a power law 
k± ~ (1.55 — t)~^'^. Hence, the small scales, in this weak turbulence formalism, are reached 
in a finite time i.e. in a catastrophic way. This is also seen on the temporal evolution of 
the enstrophy (see bottom of Figure H), with a catastrophic growth ending at t ~ 2.5, after 
which the decay of energy begins. 

Figure |12] shows the temporal evolution of the energy spectrum E±{k±, 0) of the shear- 
Alfven waves around the catastrophic time to- We see that before to, evaluated here with a 
better precision to be equal to 1.544, the energy spectrum propagates to small scales following 
a stationary k^^'^- spectrum and not a k'^- spectrum. It is only when the dissipative scale 
is reached, at to, that a remarkable effect is observed: in a very fast time the solution 
turns into the finite energy fiux spectrum /cj^ with a change of the slope propagating from 
small scales to large scales. 

Note that this picture is different from the scenario proposed by Falkovich and Sha- 
farenko (1991, hereafter FS) for the finite capacity spectra. In an example considered by FS, 
the Kolmogorov spectrum forms right behing the propagating front, whereas in our case it 
forms only after the front reaches infinite wavenumbers {i.e. dissipative region). The front 
propagation can be described in terms of self-similar solutions having a form (Falkovich and 
Shafarenko 1991; Zakharov et al. 1992): 

E^(fc^,0) = i^Eo(^), (75) 
where r = to — t. Substituting ( ffSD into the kinetic equation (|5^) we have 





which leads to the relation 

1 + 36 = a . (76) 

If Eq is stationary and has a power-law form Eq k"^, then we have another relation between 
a and b 

a + mb = (77) 
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Excluding a from ([76| ) and (^) we have l + (3+m)6 = 0. In our case this condition is satisfied 
because h = —3/2 and m = —7/3 which confirms that the front solution is of self-similar 
type. 

6 Locality of power-law spectra 

As we mentioned above, a Kolmogorov-type spectrum obtained via the Zakharov transform 
is a solution to the kinetic equations if, and only if, the original collision integral in this 
equation (before the Zakharov transformation) converges on it, - a property called locality 
of the spectrum. Having in mind that the front propagation spectrum is also of a power law 
type, let us study locality of power law spectra of a general form, 

^l(fc±,0) = A;r, EZ'{k^,0) = k'J^-% (78) 

where indices and m_s are arbitrary numbers. Recall that the collision integral in ( |54D 
is to be taken over a semi-infinite strip shown in Figure 0. It may be singular only at the 
following three points, 

(pi) K^ = L^ = +00, 

(p2) = /c_L, = 0, 
(p3) = 0, L± = k±, 

i.e. the corners and infinity of the integration area shown in Figure 1. To study conver- 
gence at point (pi) it is convenient to change variables, 

K± + L± = r+, — L± = r_, k± < < +oo, —k^ < r_ < k^. (79) 

Taking the limit +oo in the integrand (which corresponds to (pi)) and retaining the 

largest terms, we obtain the following conditions of convergence, 

m_s + rris < 0, m_s < —1- (80) 

In the vicinity of (p2) it is convenient to use the polar coordinates, 

n± = k±{l + r cos9), L± = k±rsm9, "4 ^ ^ ^ 4' -k± < r < kj_. (81) 

Considering the limit r ^ and integrating over 6 one can see that the collision integral 
converges if, and only if, > —3. Similarly, one obtains the convergence condition at point 
(p3) which is m_s > —3. 

All the convergence conditions in the kinetic equation for E"'' are, of course, symmetric 
to the case of E'^; one simply has to exchange m_s and in these conditions. Summarizing, 
one can write the conditions for simultaneous convergence for both E^ and E~'^, 

-3<m±<-l. (82) 
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The Kolmogorov spectral exponents lie on the line m+ + m„ = —4, and the locality interval 
(—3,-1) for one of them maps exactly onto the same interval for another exponent. In 
particular, the symmetric —2 Kolmogorov spectum is local. One can also see that the front 
solution with index —7/3 is local according to the above locality condition. 



7 Fokker-Planck approximation 

In the previous section, we established the fact that both the Kolmogorov —2 and front —7/3 
spectra are local. However, during the initial phase of the turbulence decay, turbulence may 
be very nonlocal. Namely, the nonlinear interaction for short waves will be dominated by 
triads that involve a long wave corresponding to the initial large-scale turbulence. Our 
locality analysis suggests that this will happen when the slope of the large-scale part of 
the spectrum is still steeper than —3, i.e. neither a —7/3 nor a —2 small-scale tail has 
grown strong enough in amplitude yet for the local interaction to take over. Further, the 
locality analysis suggests that the dominant contribution to the collision integral in this 
case will come from small vicinities of the points (p2) and (p3) both of which involve one 
small wavenumber : L and k correspondingly. Thus, one can expand the integrand of the 
collision integral in powers of these wavenumbers and reduce the kinetic equation to a second 
order Fokker-Planck equation, similarly to the way it was done for the Rossby three-wave 
process in Balk et al. (1990a, 1990b). Below, we will derive such an equation considering 
contributions of the points (p2) and (p3) separately. 

Below in this section, we will consider only the two-dimensional symmetric case and, 
therefore, we omit the superscript s in and the subscript ± for the wavevectors. The 
kinetic equation (|5^ ) can be rewritten as 



dE{k) 



dt 

where 

F(/c, K, L) = — — (cos0)^ \sin(f)\ — -— E(k 

20o 



2 J F{k,K,L)dKdL, (83) 
E{L) E{k) 



4) 



and is an angle between wavevectors k and L, so that 

cos0=^^A_Z^. (85) 

In a small vicinity of (p2) one can expand F in powers of small L and h = k — n = 0{L). 
Taking into account that cos^ (p = {^/^Y + 0{L), we have 

F{k,K,L) = '^{h/Lf\l-{h/Lf\^/^LE{k)E{L). (86) 
2oo 
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Substituting this expression into (^3|) and integrating over h from —L to L we have the 
following contribution of the point (p2), 

(87) 

where the constant D is 

Let us consider now the contribution to the collision integral that comes from the vicinity 
of the point (p3). Introducing the new variable I so that L = k + l and applying the Zakharov 
transformation (simply / —I near point (p3)), we rewrite (|83D as follows 



dE{k) 



/ dK dl[F{k,K,k + l) + F{k,K,k -I)]. (89) 

Jo J-K 



dt 

Assuming that k and I are small and that they are of the same order near point (p3) we have 



d{E{k')/k') 



+ 0{l'), (90) 

J k'=k+l/2 



E{k + l) _ E{k) 
k + l k \_ dk 

k^L^ = k^{k + lf = {k + l/2y + 0{P), (91) 

COS0 = 1 + 0{P), (92) 

I sin 01 = ^J\P - K^\{{k + 1/2)-' + 0{P)). (93) 

Substituting these expressions into (^9]) and further Taylor expanding the integrand and 
integrating over / we have the following main order contribution of the point (p3), 



where the diffusion constant D is given by (0). Combining (|8^) and (p^ we have the 
following kinetic Fokker-Planck equation, 

dE(k) d ,.d ,Eik)^^ 

The first term in the RHS of this equation conserves energy of the short waves because 
the three-wave interaction near (p3) does not transfer any energy to (or from) the large- 
scale component. The second term is not conservative: it describes a direct nonlocal energy 
transfer from long waves to the short ones. According to (p5|), the total energy of short waves 
grows exponentially. Indeed one can rewrite this equation in the form of a local conservation 
law for = e~'^^^E as follows, 

ON ^ d ,N^^ , , 
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It is interesting that this equation can also be rewritten in the form of a conservation for 
N/k\ 

^-|(-|(-)'- 

Further, there are two independent power law solutions to (pGf ) and (p7| ) : N k and 
~ 1/k. The first of these solutions corresponds to the equipartition of and a constant 
flux of N/k"^, whereas the second one corresponds to the equipartition of N/k"^ and a constant 
flux of A^. This property of the Fokker-Planck equations to have two independent integrals of 
motion such that the constant flux of one of them corresponds to the equipartition of another 
one (and vice versa) was recently noticed in Nazarenko and Laval (1998) in the context of 
the problem of passive scalars. Note however, that one could expect solutions N ~ k"^^ only 
in a very idealized situation when a short wave turbulence is generated by a source separated 
from the intense long waves by a spectral gap, and only for a limited time until the —7/3 tail 
growing from the large-scale side will fill this spectral gap. In general, the dynamics given 
by the Fokker-Planck equation (^) describes more complex combination of the instability 
and diffusion processes with an energy in-flux from the initial large scales. 



8 MHD turbulence without an external magnetic field 

We have considered until now a turbulence of Alfven waves that arises in the presence of a 
strong uniform magnetic field. Following Kraichnan (1965), one can assume that the results 
obtained for turbulence in a strong external magnetic field are applicable to MHD turbulence 
at small scales which experience the magnetic field of the large-scale component as a quasi- 
uniform external field. Furthermore, the large-scale magnetic field is much stronger than the 
one produced by the small-scales themselves because most of the MHD energy is condensed at 
large scales due to the decreasing distribution of energy among modes as the wavenumbers 
grow. In this case therefore, the small-scale dynamics consists again of a large number 
of weakly interacting Alfven waves. Using such a hypothesis and applying a dimensional 
argument, Kraichnan derived the A;~^/^ energy spectrum for MHD turbulence. However, 
Kraichnan did not take into account the local anisotropy associated with the presence of this 
external field. In Ng and Bhattacharjee (1997) (see also Goldreich and Sridhar 1997) the 
dimensional argument of Kraichnan is modified in order to take into account the anisotropic 
dependence of the characteristic time associated with Alfven waves on the wavevector by 
simply writing 

^ ~ riT • (98) 
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In that way, one obtains a kj^ energy spectrum, which agrees with the analytical and numer- 
ical results of the present paper for the spectral dependence on k±. On the other hand, the 
dependence of the spectra on k\\, as we showed before in this paper, is not universal because 
of the absence of energy transfer in the k\\ direction, although it is shown in Kinney and 
McWiUiams (1998) that for a quasi-uniform field as considered in this section, there is some 
transfer in the quasi-parallel direction. In the strictly uniform case, this spectral dependence 
is determined only by the dependence on of the driving and/or initial conditions. 

For large time, the spectrum is almost two-dimensional. The characteristic width of the 
spectrum in k\\ (described by the function fi{k\\)) is much less than its width in k±, so that 
approximately one can write 

e'{k)^CkX^5{k\i), (99) 

where C is a constant. The kj^ factor corresponds to the E± oc kj^ Kolmogorov-like 
spectrum found in this paper (the physical dimensions of and being different) . In the 
context of MHD turbulence, this spectrum is vahd only locally, that is for distances smaller 
than the length-scale of the magnetic field associated with the energy-containing part of 
MHD turbulence. Let us average this spectrum over the large energy containing scales, that 
is over all possible directions of Bq. Writing k± — |k x Bo|/|-Bo| and A;|| = |k • 5o|/|-Bo| and 
assuming that Bq takes all possible directions in 3D space with equal probability, we have 
for the averaged spectrum 

(e^ = I e^(k, x) da{C) = J CS{C • k) |C x k\-'da{C), (100) 

where C = (sin 9 cos 0, sin 9 sin (p, cos 9) is a unit vector in the coordinate space and da — 
sin 9 d9 dcf) is the surface element on the unit sphere. Choosing 9 to be the angle between k 
and Bq and (f) to be the angle in the transverse to k plane, we have 

{e^)^Crdct> r^-^^^^su.9d9^2nCk-\ (101) 

Jo Jo \k\ \k\'^ 

This isotropic spectrum represents the averaged energy density in 3D wavevector space. By 
averaging over all possible directions of the wavevector, we get the following density of the 
energy distribution over the absolute value of the wavevector, 

Ek^Sn^Ck-^. (102) 

As wc see, taking into account the local anisotropy and subsequent averaging over the 
isotropic energy containing scales results in an isotropic energy spectrum k~'^. This re- 
sult is different from the spectrum derived by Kraichnan without taking into account 
the local anisotropy of small scales. The difference in spectral indices may also arise from 
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the fact that the approach here is that of weak turbulence, whereas in the strong turbulence 
case, isotropy is recovered on average and a different spectrum - that of Kraichnan - obtains. 

Solar wind data (Matthaeus and Goldstein 1982) indicates that the isotropic spectrum 
scales as with a ~ 1.67, close to the Kolmogorov value for neutral fluids (without 
intermittency corrections which are known to occur) ; hence, it could be interpreted as well 
as being a Kraichnan-like spectrum steepened by intermittency effects which are known to 
take place in strong MHD turbulence as well in the form of current and vorticity filaments, 
ribbons and sheets, and magnetic flux tubes. However in the context of the interstellar 
medium (ISM), data show that the velocity dispersion is correlated with the size of the 
region observed (Larson 1981; Scalo 1984; Falgarone et al. 1992). These correlations are 
approximately of power-law form such that the corresponding energy spectrum scales with 
a spectral index ranging from a ~ 1.6 to a ~ 2. Then the weak turbulence approach could 
explain the steepening of the spectrum. But the variety of physical processes in the ISM, 
such as shocks or dispersive effects for instance, do not allow to give a definitive answer but 
rather to ask the question : by how much is the energy spectrum of a turbulent medium 
affected by such physical processes ? A formalism that incorporates dispersive effects in 
MHD, e.g. the Hall current for a strongly ionized plasma, as in the magnctosphcrc or in 
the vicinity of proto-stellar jets or the ambipolar drift in the weakly ionized plasma of the 
interstellar medium at large, will be useful but is left for future work. So is the incorporation 
of compressibility. 

9 Conclusion 

We have obtained in this paper the kinetic equations for weak Alfvenic turbulence in the 
presence of correlations between the velocity and the magnetic field, and taking into account 
the non-mirror invariance of the MHD equations leading to non-zero helical terms. These 
equations, contrary to what is stated in Sridhar and Goldreich (1994), obtain at the level of 
three-wave interactions (see also Montgomery and Matthaeus 1995; Ng and Bhattacharjee 
1996). 

In this anisotropic medium, a new spectral tensor must be taken into account in the 
formalism when compared to the isotropic case (which can include terms proportional to the 
helicity) ; this new spectral tensor P is linked to the anisotropy induced by the presence of a 
strong uniform magnetic field, and we can also study its dynamics. This purely anisotropic 
correlator was also analysed in the case of neutral fluids in the presence of rotation (Cambon 
and Jacquin 1989). 

We obtain an asymptotic two-dimensionalisation of the spectra: indeed, the evolution 
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of the turbulent spectra at each k\\ is determined only by the spectra at the same /cy and 
by the purely 2D state characterized by k\\ = 0. This property of bi-dimensionalisation 
was previously obtained theoretically from an analysis of the linearised MHD equations 
(Montgomery and Turner 1981) and using phenomenological models (Ng and Bhattacharjee 
1996), and numerically as well (Oughton et al. 1994; Kinney and McWilliams 1998), whereas 
it is obtained in our paper from the rigorously derived kinetic equations. Note that the strong 
field Bq has no structure (it is a k = field), whereas the analysis performed in Kinney 
and McWilliams (1998) considers a strong quasi-uniform magnetic field of characteristic 
wavenumber fcj^ 7^ 0, in which case the authors find that bi-dimensionalisation obtains as 
well for large enough wavenumbers. 

We have considered three-dimensional turbulence in the asymptotic regime of large time 
when the spectrum tends to a quasi-2D form. This is the same regime for which RMHD 
approach is valid (Strauss 1976). However, in addition to the shear- Alfven waves described 
by the RMHD equations, the kinetic equations describe also the dynamics of the so-called 
pseudo-Alfven waves which are decoupled from the shear-Alfven waves in this case, from 
the magnetic helicity and from the pseudo-helicity. Finding the Kolmogorov solution for the 
3D case is technically very similar to the case of 2D turbulence. This leads to the spectrum 
/(/c||) fcj^, where f{k\\) is an arbitrary function which is to be fixed by matching in the forcing 
region at small wavenumbers. The k\\ dependence is non-universal and depends on the form 
of the forcing because of the property that there is no energy transfer between different fcy's. 

The f{k\\) k~^ energy spectrum is well verified by numerics, which also shows that this 
spectrum is reached in a singular fashion with small scales developing in a finite time. We 
also obtain a family of Kolmogorov solutions with different values of spectra for different wave 
polarities and we show that the sum of the spectral exponents of these spectra is equal to 
—4. The dynamics of both the shear-Alfven waves and the pseudo-Alfven waves is obtained. 
Finally, the small-scale spectrum of isotropic 3D MHD turbulence in the case when there is 
no external field is also derived. 

The weak turbulence regime remains valid as long as the Alfven characteristic time 
(fc||6o)~^ remains small compared to the transfer time (which can be evaluated from equation 
(P^)), that is to say for E']_{k^,Q)k\/ Bl <^ k\\/k^. Using the exact scaling law found in 
this paper, £'^(/c^,0) oc A;J^, we see that the condition for weak turbulence is less satisfied 
for large /c^, or small k\\. This means that we have a non-uniform expansion in Bq. 

The dynamo problem in the present formalism reduces to its simplest expression : in the 
presence of a strong uniform magnetic field Bq, to a first approximation (closing the equations 
at the level of second-order correlation tensors), one obtains immediate equipartition between 
the kinetic and the magnetic wave energies, corresponding to an instantaneous efficiency of 
the dynamo. Of course, one may ask about the origin of Bq itself, in which case one may 
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resort to standard dynamo theories (see e.g. Parker 1994). We see no tendency towards 
condensation. 

In view of the ubiquity of turbulent conducting flows embedded in strong quasi-uniform 
magnetic fields, the present derivation should be of some use when studying the dynamics 
of such media, even though compressibility effects have been ignored. It can be argued 
(Goldreich and Sridhar 1995) that this incompressible approximation may be sufficient be- 
cause of the damping of the fast magnetosonic wave by plasma kinetic effects. Note however 
that Bhattacharjee, Ng and Spangler (1998) found that in the presence of spatial inhomo- 
geneities, there are significant departures from incompressibility at the leading order of an 
asymptotic theory which assumes that the Mach number of the turbulence is small. Finally, 
the wave energy may not remain negligible for all times, in which case resort to phenomeno- 
logical models for strong MHD turbulence is required. Is desirable as well an exploration of 
such complex flows through analysis of laboratory and numerical experiments, and through 
detailed observations like those stemming from satellite data for the solar wind, from the 
THEMIS instrument for the Sun looking at the small-scale magnetic structures of the pho- 
tosphere, and the planned large array instrumentation (LSA-ALMA) to observe in detail 
the interstellar medium. 
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Appendix 



From the dynamical equations (^4[) one writes successively for the second and third-order 
moments of the fields : 

-tek'^PM^') J {a-'(/c)a^'(L)a|(k)} e-2-'-*5«L,k'4L 

and 

dt{a^^{k)a^:{k')af{k")}= (104) 

-zek^P,n{k) J {a„^(«)<(L)a^<(k')a^^;;(k")}e-2*^-'=*<5„L,k4L 

+ {[k,s,j] ^ [k',s',f] ^ [k",s",f] ^ [k,s,j]} . 

Asymptotic closure for the leading order contributions to each of the cumulants follows 
from the following procedure or algorithm. Cumulants are in (1:1) correspondence with the 
moments : in the zero mean case, the second and third moments are the second and third 
cumulants; the fourth order moment is the sum of cumulants where each decomposition of the 
moment is taken once, namely the sum of the fourth order cumulant plus products of second 
order ones. One attempts to solve the hierarchy of cumulant equations perturbatively. The 
asymptotic expansions generated in this way are not uniform in time because of the presence 
of small divisors (mainly but not totally due to resonances). In order to restore the uniform 
validity of the asymptotic expansions we must allow the leading order contributions to each 
of the cumulants vary slowly in time and choose that time dependence to achieve that goal. 
Where necessary, and where the notion of well-orderedness does not make sense in Fourier 
space, one must look at the corresponding asymptotic expansions in physical space. The 
result is another set of asymptotic expansions which include both the kinetic equations for 
the second order moments, combinations of which give the parallel and total energies and 
helicity densities, and similar equations for higher order cumulants which can be collectively 
solved by a common frequency renormalization. The kinetic equations are valid for time 
scales of the order of and possibly longer depending on how degenerate the resonant 
manifold are. Success in obtaining asymptotic closure depends on two ingredients. The first 
is the degree to which the linear waves interact to randomize phases and the second is the 
fact that the nonlinear regeneration of the third order moment by the fourth order moment 
in equation (|104| ) depends more on the product of the second order moments than it does 
on the fourth order cumulant. We now give details of the calculations. 

The fourth-order moment in the above equation, (/cLk'k") in short-hand notation, de- 
composes into the sum of three products of second-order moments, and a fourth-order 
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cumulant {mnj'j"}. The latter does not contribute to secular behavior, and of the remain- 
ing terms one is absent as well in the kinetic equations because it involves the combination 
of wavenumbers (kL) (k'k") : it introduces, because of homogeneity, a 6{k + L) factor which 
combined with the convolution integral leads to a zero contribution for k = 0. Hence, the 
time evolution of the third-order cumulants leads to six terms that read : 



t\„2isui"t 



2is'Lu"t 



d, {a^^{k)a^:{k')ap{k")} = -zefc„P,„(k)g„^.?'(k')g:5':(k")e2--'* 
-^efc^P,„(k)g-^f;' (k")g^5:(k')e 

-zeCP,"n(k")g-^^'(k')g£^(k)e 



(105) 



2is"Lu't 



It can be shown that, of these six terms, only the fourth and fifth ones give non-zero 
contributions to the kinetic equations. Defining 



and integrating equation ( p.05| ) over time, the exponential terms will lead to 

exp [iuJk,KL\ - 1 



exp [ituJk,nL]dt 



(106) 



(107) 



Substituting these expressions in (|105|) , only the terms which have an argument in the A 
functions that cancel exactly with the arguments in the exponential appearing in (|103|) will 
contribute. We then obtain the fundamental kinetic equations for the energy tensor, viz. : 



dtqf;,{k')5{k + k') 



(108) 



-e^kmPjnik) J A:2pP„g(L)gp^-^(K)gg^j,'(k')A(-2sc^i)5«L,kC?KL 
-e^krnPW^)S,s' J A;;P,,,(k')g;^'-^(K)g,l^(L)A(-2su;i)5«L,k4L 
-e'kiP,>^{k') J k,,P^,{L)q;::-^\K)q;f{k)A{-2s'u;,)S,j^yd,^ 

-e^k'^Pj,nik')Sss' I A;pP„g(k)g-^-^'(/c)9g"'(L)A(-2scui)5«L,k'4L • 
We now perform an integration over the delta and taking the limit t +oo we find 

dt[qg,ik')6ik + k')]= (109) 

/ /{Qr(«)^.n(k)p,.(L)[g-:(k') 

+gr'(«)^/n(k)p„KL)k^'^(k) 



c-/ N 1 
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-7r5,,,g,^(K)P,v,(k)P,„(k')g,::(L)5(/s:||)|t/«:irf/€2t^/t||, 

where V stands for the principal value of the integral. 

In the case where s = s' of interest here because the cross-correlators between z-fields of 
opposite polarities decay to zero in that approximation, the above equations simplify to : 

f9ag-,(k')±gf_,.(k')]= (110) 

2e2 / P,„(k)P,,,(k)[g-(L)±g-(L)]Q-(«)5(«:|,)rf«:irf«:2^^«:|| 

P,„(k)P„,(L)[g^^j,(k') ± q^f^{y)]Q-^\K)5{K^)dK^dK2dK\\ 

-e^J P,v„(k)P„,(L)[g|^^(k') ±g,^j(k')]Qr(«)<^('«||)rf'^irf'«2C^«:|| 

+ / P,„(k)P„,(L)[g-,(kO T q]W]Q,^iK /''''^''''^''^^ 

SIX J /t^II 

-^p / p,..(k)p„,(L)[g-(ko Tg-(kO]gr(/^) '^"^t''^"'' • 

SIX J rb|j 

To derive the kinetic equations we need now to develop 

aie^(k) = dMl{{^) + gi(k) + gf|(k)) , (111) 



ai$^(k) = kl^9igp(k) 
1 

-2^k^ 

~4; 



9iP(k) = \^i%{kMn{^) + ?p(k)) - h{q^\{\^) + gp(k))) , 

in terms of the above expressions ( |109| ) and ( |110| ). This leads to : 

dte'{\^) = (112) 
^ j [2(g-(L) + g-(L) + gf|(L) - g-(k) - g-(k) + ^^(k)) 

+ ^(C(k) + g^Kk)) - ^(C(L) + Qnlim Q,'if^)SiK\\)dK^dK2dK\\ , 
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9t$^(k) = (113) 
^ / [2gf|(L) - 2*(g|(L) + 5ff(L)) + Sfc„A;KC(L) + ?^KL))] 
+[-2g||(k) + (Ml + Mi _ ^fci^A)(g|(k) + gff(k))] Q-(«)(^(«||)d«:id«:,d«:|| 



ati?^(k) = (114) 
|^[2 / (g-(L) - ?^2(L)) - ^(C(L) - 9if (L)) - ^(C(L) " C2(L)) 

+ ^I^(C(L) - gnKL)) Ql'{K)5{K\{)dKidK2dK\\ 

- J {Qtm - - ^(c(k) - ?ir(k)) - ^{Qtm - 

+ ^^^^iQtm - ?lf(k)) Q-,'{K)S{K\\)dK,dK2dK\\ 

- J (Qtm - Qiim - ^{Qim - Qtim - ^m^) - c(k)) 

+ ^^^^^(lni^) - Ini^)) Qk"{'^)H'^\\)dKidK2dKii 
/ [-T^(^^^2^(k) + g^I(k)) - ^(gfl(k) + g^f (k)) 
+^^^^(?r2 (k) + (k)) + ^(g^f (k) + g-(k)) + ^(g^f (k) + C(k)) 

'^^•'^^-^(..7(k) + .,T(k))]gr(«)"^^''^^''^" 



k2L2 
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9tP(k) = (115) 

2fc1 



.[2 1 (g,^|(L) + g^^(L)) - ^(g|(L) + g^f(L)) - -L!^(gi-(L) + g-(L)) 



■ (C(L) +g:KL))Qfe'(K)5(/t|i)c?/tirf/t2c?«;ii 

?(k) + g^|(k)) + (-^1^ - ^ + (k) + gi7(k)) (K)5('^||)rf'^irf'^2rfK|| 

(g^^(k) + g^|(k)) + (-:^ - ^ + (k) + g|(k)) g,- («)'5('^||)^'^irf'^2rf«:|| 

isP f, L\\Lj knh fciik • LL 



+(-^ - ^ + ^^^)(C(k) - .i(k)) Qr(«)'^^"^^'^" 



^[2 / (gp(L) + g^^KL)) - ^(gf^^lL) + g-(L)) - ^{q^^) + g^KL)) 



k2k\\knki 



k4 



L2 k2 k2L2 - K\\ 

The final step which leads to the kinetic equations (|26|)-(p9D, consists in substituting the 
expressions (^) inside ( |112| )-( PT5D . For this computation it is useful to note that : 

X2 + F2 = ^^^2 ^ (IIQ) 

X2 + = klLl , 
Z^-X^ = {Ll - Ll){kl - kl) + Akik2L^L2 , 
XZ = LiL2{kl - kl) + kik2{Ll - Ll) , 
XY = kl{k2Li - kiL2) + LiL2{kl - kl) + kik2{Ll - L\) . 
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Figure 1: Geometrical representation of the Zakharov transformation. The rectangular 
region, corresponding to the triad interaction = L + /c, is decomposed into four different 
regions called (1), (2), (3) and (4); uji and uj2 are respectively the dimensionless variables 
K±/k± and L±/k±. The Zakharov transformation applied to the collision integral consists in 
exchanging regions (1) and (2), and regions (3) and (4). 
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Figure 2: Variation of \JCK{ns) CK^—ng — 4) as a function of — n^. Notice the symmetry 
around the value —Ug corresponding to the case of zero velocity-magnetic field correlation. 
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Figure 3: Variation of P±/P± , the ratio of fluxes of energy, as a function of —Ug. For the 
zero cross correlation case the ratio is 1. 
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Figure 4: Variation of ^JC'j^{ns) C'j^{—ns — 4) as a function of —Ug- Notice the symmetry 
around the value —Ug — 2 corresponding to the zero cross correlation case. 
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Figure 5: Temporal evolution of the energy E±(t) (top) and the enstrophy {u}^[t)) in units 
of 1. 10^. Notice the conservation of the energy up to the time to ~ 1.55. 
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Figure 6: Energy spectra E±{k±,0) of the shear-Alfven waves in the zero cross correlation 
case for the times t — (dot), t — 1.0 (dash-dot), t — 1.5 (small-dash), t — 1.6 (sohd) and 
t — 10.0 (long-dash) ; the straight line follows a kj^. 
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Figure 7: Temporal evolution of energies (top) E'l (solid) , £'J (long-dash) , E'^ (small-dash) 
and (dash-dot) ; the same notation is used for the enstrophies (bottom) which are in 
units of 1. 10^. Notice that energies are conserved till the time ~ 5. 
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Figure 8: Temporal evolution of the cross correlations p± (solid) and py (dash). These 
quantities are conserved up to the time ^q. 
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Figure 9: Compensated spectra E]^E]^kj_ at times t — (sohd), t — A (long-dash), t — 6 
(small-dash), t — 8 (dash-dot) and t — 20 (dot). 
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Figure 10: Compensated spectra E^E^k\ for the same times and with the same symbols as 
in Figure |. 
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Figure 11: Temporal evolution (top), in lin-log coordinates, of the front of energy propagating 
to small scales. The sohd line and the dash-dot line correspond respectively to an energy of 
10~^^ and 10~^^. An abrupt change is visible at time to ~ 1.55 (vertical dotted hne). The 
log(A;_L) as a function of log(1.55 — t) (bottom) displays a power law in k± ~ (1.55 — t)~^-^ 
(large dash hne). 
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Figure 12: Temporal evolution of the energy spectrum E_i_{k±,0) of the shear-Alfven waves 
around the catastrophic time to ~ 1.544. For t < to (top), t = 1.50 (dot), t = 1.53 (dash- 
dot), t — 1.54 (dash), t — 1.542 (long dash) and t — 1.543 (solid)) a k]^'^^^- spectrum is 
observed. For t > to (bottom), i = 1.544 (solid), t = 1.546 (long dash), t = 1.548 (dash), 
t — 1.55 (dash-dot) and t — 1.58 (dot)) a fast change of the slope appears to give finally a 
/cj^- spectrum. Note that this change propagates from small scales to large scales. In both 
cases straight lines follow either a kj^^^ or a kj^. 
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